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Abstract 

SICANE is a neutron scattering multidetector facility for the determina- 
tion of the quenching factor (ratio of the response to nuclear recoils and to 
electrons) of cryogenic detectors used in direct WIMP searches. Well colli- 
mated monoenergetic neutron beams are obtained with inverse (p, n) reac- 
tions. The facility is described, and results obtained for the quenching factors 
of scintillation in Nal(Tl) and of heat and ionization in Ge are presented. 



1 Introduction 



Since 1985 [|IJ, various experimental techniques have been proposed to detect 
the nuclear recoils that would be produced by the scattering of WIMP dark 
matter particles from our galactic halo (for a recent review of experimental 
searches, see e.g. Ref. 0). The goal is to be sensitive to recoil energies 
down to 10 keV ||, at event rates well below one per kilogram of detector 
material per day For this, the detector response to low-energy nuclear 
recoils must be well known. Depending on the detection process (ionization, 
scintillation or heat), this response may differ significantly from that inferred 
from calibration with electron or gamma-ray sources. A quantity of interest 
is thus the quenching factor, defined as the ratio of the signal amplitudes 
induced by a nuclear recoil and an electron of the same energy. This factor 
depends mainly on the detector material, the energy of the recoil and the 
detection process, although temperature and alignment effects may play a 
role (see e.g. Ref. 0). 

Quenching factor measurements have been made on scintillation detectors 
such as Nal ||, Csl 0, CaF 2 [§] and Xe ||, leading to values ranging from 
4 to 30% depending on the recoil nucleus and its energy. The measured 
quenching factors of ionization detectors such as Si |TU] or Ge [Tl"| , |T2| 



vary between 25 and 40% and appear to follow the recoil energy dependence 
predicted by the Lindhart theory |16[ . More recently, the thermal detection 



efficiency for recoiling nuclei has been measured for the first time in a TeO"2 



cryogenic detector, giving a thermal quenching factor slightly above unity []TT 



The purpose of the SICANE facility (Site de CAlibration NEutron) is the 
measurement of quenching factors of cryogenic detectors using nuclear recoils 
induced by monoenergetic neutron beams. The facility has been commis- 
sioned with tests using scintillation in a Nal(Tl) crystal and measurements 
of quenching factors for ionization and heat signals in a germanium cryogenic 
detector developed by the EDELWEISS collaboration. 
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2 Experimental Procedure 



2.1 Principle and Motivation 



The basic method for the measurement of the quenching factor of a detector is 
the following. First the detector response to electron recoils is obtained using 
calibrated 7-ray source. This yields the calibration of the signal amplitude 
in terms of an electron- equivalent energy, E ee . In a second step, the detector 
is exposed to a source producing nuclear recoils with a known kinetic energy 
Er. The quenching factor Q is then 

Q = ^ (i) 

The main difficulty for a precise and efficient measurement is the production 
of recoils with known En values. 

In the present method, nuclear recoils are produced by neutron scattering. 
Other methods for producing nuclear recoils exist, such as the implantation 



of a-particle sources close to the surface of the detector [17], [Tj|. When the 
a escapes the detector, the response of the detector to the recoiling daughter 
can be measured accurately. However this method is not adequate for detec- 
tors which response differs according to whether the interaction occurs close 



to surface or in the bulk of the volume fll9 |. Another drawback is that the 
recoiling daughter is often very different from the detector material. Neutron 
interactions have the advantage of occurring rather uniformly throughout the 
detector volume. The recoil nuclei thus produced are very similar to those 
that would be induced by WIMP scattering. 

In Ge cryogenic detectors with simultaneous read-out of the heat and 
ionization signals the recoil energy can be measured event-by-event, inde- 
pendently of the nature of the recoil (electron or nuclear). In this type of 
detector, it is possible to directly measure the ratio of the ionization and heat 
quenching values[| as a function of recoil energy |13] . If one assumes that the 



heat quenching factor Q' is unityn, these data can be interpreted as ionization 



■'■For clarity, from now on the symbol for ionization or scintillation quenching factors 
will be Q, and Q' will be used to represent heat quenching factors. 

2 This assumption is reasonable if the detector can be approximated as a close system 
where all the excitation energy is thermalized within a time scale inferior to that of the 
read-out of the heat signal and no energy is stored in crystal defects or other long-life 
processes. 



2 



quenching measurements. Alternatively, these data can be compared to the 
available direct measurements of ionization quenching in order to provide a 
verification of how close to unity are the heat quenching factors ||13|| . This 
comparison assumes that the ionization quenching factors do not depend on 



temperature, as all direct measurements |L1], [12], |14[ were performed at liquid 



nitrogen temperature (77 K) instead of cryogenic temperatures (~25 mK). 

This method cannot be applied to detectors with a single read-out, or to 
heat-and-scintillation detectors, as the amount of energy escaping as light is 
difficult to control. In addition, scintillation yields are known to depend on 
crystal temperature (see e.g. Ref. |0| ). 

Clearly, a method providing nuclear recoils of a known energy is relevant 
for the calibration of cryogenic detectors. In a non-relativistic elastic neutron 
scattering, the energy of the recoil nuclei Er is given by 

F (ff\ - F 4Mm n (1-COS0) 

Er{0) ~ En {m n + MY 2 (2) 

where E n is the incident neutron kinetic energy, 9 is the neutron scattering 
angle in the center-of-mass frame and m n and M are the neutron and nucleus 
masses, respectively. Even with a monoenergetic neutron beam, Er values 
vary from zero to _Er(180°). In experiments where 9 is not measured some 
information about the quenching factor can be obtained from the maximum 
E ee value observed for a given E n , the quenching factor being then equal 
to E ee / Er(180°) [11]. A more precise measurement consists in measuring 



the scattering angle of the neutron JTI]. This requires neutron detectors 



with a small solid angle. However, in order to be efficient the total solid 
angle for the detection of the scattered neutrons must be as large as possible. 
This is a crucial point for performing in a reasonable time quenching factor 
measurements of cryogenic heat detectors, as their response is intrinsically 
slow. 



These precision and efficiency constraints naturally lead to a concept of a 
large array of neutron detectors. With 48 neutron detectors, SICANE is the 
largest existing array. With such an extensive array, scattering at different 
angles can be measured simultaneously, reducing the systematic errors, and 
improving the diagnostic tools to identify the various background sources. 

The basic performance of this array can readily be tested using scatter- 
ing on a small Nal detector at room temperature, for which scintillation 
quenching data already exist 0. For actual measurements on cryogenic de- 
tectors, two new potential problems arise. Firstly, the cryostat surrounding 
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the detector represents an important amount of additional material where 
neutrons may scatter before or after interacting in the detector itself. Sec- 
ondly, cryogenic detectors are intrinsically slow, and some care must be taken 
in establishing clean coincidences with the relatively faster array of neutron 
detectors while retaining a sufficiently high count rate. For the commission- 
ing measurements of an actual cryogenic detector, a Ge heat-and-ionization 
detector was chosen as the experimental situation for these detectors is rather 
well controlled (see above), while still awaiting direct measurements of the 
low temperature behavior of both Q and Q'. 



2.2 Neutron Beams 



Up to now, experiments involving these techniques used monoenergetic neu- 
trons produced by deuteron beam interactions on a deuterated target, or by 
a proton beam interacting on a 7 Li target fll4 |. The neutron emission in 



these reactions is almost isotropic. This reduces the flux impinging on the 
target detector, and imposes that the scattered-neutron detector array be 
very carefully shielded from direct neutrons [2D|, putting severe constraints 
on the spatial extension of the array. 

In contrast, reactions with inverse kinematics (heavy ion beam on hy- 



drogenated target) produce a naturally collimated neutron emission [^T], |22 
Using endothermal reaction enhances the focusing of the neutron emission in 
the forward direction. The neutron flux on the detector under study is thus 
increased while the flux on the surrounding array is reduced. The neutrons 
are all produced within a small opening angle Q max given by: 

8in2e _ = ^(l-^) (3) 
m n Mi Ei 

where is the threshold energy, Ej is the incident kinetic energy, and Mj, 
m p , m n and Mf are the masses of the projectile, the proton, the neutron 
and the residual nucleus, respectively. In terms of the (negative) Q-value of 
the reaction, E t h is equal to —Q(Mi + m p )/m p . At a given angle relative to 
the beam axis < Q m ax, there are two monoenergetic neutron groups with 
energies E n \ and E n2 given by: 

2QMi 



E nl , n2 = , ~~~T""T7\ {x + cos 2 6 ± Jcos 2 Q{x + COS 2 9)} (4) 

x[m v + Mi) v 

where x = E t h / Ei, and terms of order Q/Mi and (m n — m p )/Mj have been 
neglected. 
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The optimization of the neutron beam has been performed in experiments 
conducted at the tandem accelerator facilities of the Institut de Physique 
Nucleaire d'Orsay. Beams of 7 Li [^TJ and n B [B^] have been tested using an 
oscillating hydrogenated polyethylene target with a thickness of 200 /ig-cm -2 . 
The beams were pulsed at 2.5 MHz with 5 ns wide bursts. 

The contents of the induced neutron beams were studied by measuring 
the time-of-ffight between the beam bursts and the arrival in a NE213 scin- 
tillation detector located at an angle of O varying from to 25 degrees at 
a distance of 110 cm away from the hydrogenated target. The time-of-flight 
spectra recorded at a 11 B beam energy of 35 MeV for different values of 
are shown in Fig. [l|. The two neutron groups (labeled n\ and 712) and the 
prompt 7 produced in reactions in the hydrogenated target (labeled 7) are 
clearly seen. The positions of the n\ and n<i peaks follow the expectation of 
Eq. f|, as illustrated in Fig. |2|, merging together when the maximum value 
©max = 13.1° is reached. In addition, an important background under the 
rii and n 2 peaks is observed in Fig. [I]. This continuum remains even beyond 
the Q m ax value imposed by the kinematics of the inverse (p,n) reaction. It 
is also present at all beam energies. It is therefore attributed to neutrons 
induced in fusion reactions on 12 C nuclei from the polyethylene target. 

As shown in Fig. [| this background is found to be less important with the 
7 Li beams, which were therefore selected for the following experiments. The 
selected 7 Li beam energies for the Nal and Ge measurements are 13.7 and 14 
MeV, corresponding to n 2 neutron kinetic energies at O=0° of 2.15 and 2.32 
MeV and maximum emission angles Q m ax of 11.5° and 13.4°, respectively. 

The neutron beam energy resolution is dominated by the energy losses 
in the target (150 keV over the entire target thickness). The contributions 
from the Li beam properties themselves and straggling are 14 and 4 keV, 
respectively. The resulting uncertainty on E n is 59 keV. The calculated 
neutron flux |^TJ is approximately of 10 5 n/s/str for a 7 Li beam intensity of 
4 nA. 



2.3 Neutron Detector Array 

In the quenching measurements, the NE213 detector at 6 = 0° is removed, 
and the detector under study is installed 60 cm away from the polyethylene 
target. The scattered neutrons are detected in an array of 48 NE213 scintilla- 
tors located at a distance of 100 cm from the central detector and distributed 
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in four rings corresponding to scattering angles of 45°, 90°, 120° and 165°. 
Thus, in one experiment the quenching factor is measured at four different 
values of 6, and consequently at four different values of Er. The recoil energy 
values probed in the Nal and Ge experiments are listed in Tables [I] and [| 

The NE213 hexagonal cells have an external diameter of 10 cm and a 
thickness of 4 cm. A photomultiplier is optically coupled to each cell via a 
light guide. The support structure of the 48 neutron detectors is shown in 
Fig. It is made of dural and steel, thin enough to minimize additional 
neutron scattering. The detector positioning is performed in the following 
way. First, a laser beam is propagated along the accelerator beam axis, 
simulating the neutron beam at 9 = 0°. A rotating mirror is then set at 
the center of the array. The mirror is then tilted to the desired scattering 
angle and is rotated around the beam axis. The alignment of each NE213 
detector is then adjusted so that the center of its front face is illuminated by 
the reflected laser beam. 



2.4 Electronics 

The background due to 7 interactions in the NE213 cells is reduced by using 
the neutron-7 discrimination capabilities of these detectors. The signal com- 
ing from each detector is sent to two charge digitizers, integrating the current 
during Ti=40 ns and T 2 = 400 ns, respectively. The ratio of the correspond- 
ing integrated charges q\ and g 2 provides a good 7/neutron discrimination. 

The energy signal E ee from the center detector is obtained with standard 
spectroscopy amplifiers in the cases of scintillation in Nal and ionization in 
Ge. The much slower heat signals of the Ge cryogenic detector are recorded 
with a wave form digitizer. 

For each event, the following parameters are stored: the q\ and g 2 values 
for all NE213 cells, as well as their time-of-flight relative to the beam burst, 
and the energy and time-of-flight measurements of the Ge or Nal(Tl) detector 
under study. As it will be seen, the time delay between the beam burst and 
the hit in the neutron detector is a powerful tool to identify true neutron 
scattering events. This delay Ati_3 is the sum of those incurred as the 
neutron goes i) from the polyethylene target to the center detector (Ati_ 2 ) 
and ii) from that detector to a NE213 cell (At 2 ~3)- With a fast center 
detector such as Nal, additional reduction of random coincidences may be 
achieved by a further discrimination on either Ati_ 2 or At 2 _ 3 
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An event is denned as a coincidence between the center detector and at 
least one of the 48 neutron detectors. Data read-out is triggered by coinci- 
dences between any NE213 detector and the center detector within a time 
window of 1 /xs. This width is chosen so as to accommodate the time resolu- 
tion on the lowest energy signals in the center detectors. 

The most delicate aspect of the trigger is to reconcile the speed of the 
NE213 signal read-out with the slowness of the coincidence decision. The 
solution adopted here is to initiate the digitization of the NE213 information 
as soon as any of them pass a certain threshold, and to hold this information 
until the coincidence decision with the slower Ge detector is issued. If no 
coincidence arrives, the NE213 electronics is reset. Otherwise, the NE213 
data is further held until the end of the Ge heat and ionization read out, 250 
ms later. 

For the tests with the faster Nal detector, the width of the coincidence 
with the neutron array can be reduced to 500 ns, and the stand-alone NE213 
pre-trigger is not necessary. 

2.5 Recoil Energy Resolution 

The accuracy of the determination of Er on an event-by-event basis depends 
on the uncertainty on E n and on 9. In natural Ge, an additional 2.6% rms 
fluctuation on Er comes from the difference between the isotope masses. 

The energy loss of the 7 Li ions in the hydrogenated target leads to a 
dependence of the neutron kinetic energy on the position of the reaction along 
the target length. For 2.32 MeV neutrons, the energy spread on E n is ±3%, 
leading to an equivalent spread on Er. The main cause of event-by-event 
fluctuations of 9 is the solid angle of the NE213 cells. Their contributions to 
the r.m.s. width of the E R distribution decreases from 2.5% at 90° to 0.6% 
at 165°. All these effects taken together contribute to a 4 % r.m.s. event- 
by-event fluctuation of Er. This event-by-event uncertainty is much smaller 
than the resolution that can be achieved by a time-of-flight measurement 
over a distance of 1 m. It is also significantly smaller than typical resolutions 
for the measurement of E ee in Nal(Tl) or Ge detectors. 
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3 Monte Carlo Simulations 



The experimental precision on Q arises from the possibility to select quasi- 
monoenergetic neutrons and from the simple kinematics of elastic collisions. 
The determination of Er depends only on the kinematics of the neutron- 
producing reaction (yielding E n ) and the accurate positioning of the NE213 
cells (yielding 6). 

To preserve this simplicity, it is important to reduce as much as possible 
the multiple scattering of neutrons and all sources of beam-correlated signals 
between the target detector and the NE213 cells. To investigate this, the 
experiment for the measurement of the heat and ionization quenching fac- 



tors in Ge has been simulated using the Monte Carlo program GEANT |]23 
In these simulations neutrons are produced one by one, thus assuming that 
the beam current and bunch structure is such that the pile-up of interac- 
tions coming from two different neutrons is negligible. In addition to the 
Ge detector, the simulation takes into account interactions occuring in the 
material composing 48 NE213 detectors and photomultipliers, the cryostat 
and the reaction chamber. It also takes into account of the thickness of the 
hydrogenated target and follows the flight of the neutron from its creation 
to its detection in one of the NE213 cells. A cut on the time-of-flight similar 
to the one later applied to the data is also applied to the simulated data. 

These simulations take into account the inefficiencies of the NE213 cells 
due to the experimental energy thresholds and their effect on the detection 
of the elasticaly and inelasticaly scattered neutrons (energies of 2.2 and 1.6 
MeV, respectively: see below). 

Predicted rates for elastic and inelastic scattering events, including multiple- 
scattering events, were derived from these simulations, assuming an inte- 
grated flux of 70 neutrons impinging on the Ge detector per second, corre- 
sponding to a 7 Li beam current of 5 nA. The single rate in the center Ge 
detector is approximately 10 neutrons per second. 

Table [| lists the Ge-NE213 coincident rates for the different event types 
described in the following. 

Inelastic Neutron Scattering 

In addition to elastic scattering, neutrons can have inelastic interaction 
with Ge nuclei. The end products of such a collision are i) a nuclear recoil 
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with a kinetic energy inferior to that associated with elastic scattering, ii) 
a scattered neutron with a reduced energy and thus an increased time-of- 
flight to the NE213 cells, and Hi) a 7 emission which may deposit some of 
its energy in the detector. Given the neutron energy, the most important 
7's to be expected |24]] are the 596 keV transition in 74 Ge, 835 keV in 72 Ge 
and 563 keV in 76 Ge. The simulation predicts that 82% of these 7 escape 
detection, leaving only the neutron time-of-flight as a mean to discriminate 
these events from elastic scatters. The cross-sections for each of these states 
are of the order of 1 b [24|, compared with an elastic cross-section of the 



order of 2 b EH]. Given the isotopic abundances in Ge, the 563 keV yield is 



only a fifth of the 596 keV one. 

The recoil energies produced by these inelastic processes are approxi- 
mately 6 keV inferior to the values for elastic scattering. The differen- 
tial cross-section of the inelastic process has a pronounced maximum for 
backward-scattered neutrons, and it is only at these angles that it can result 
in an appreciable background. At 120°, the coincident rate of inelastic events 
is approximately half of the elastic rate (Table ||]). These two types of events 
can be identified by their ~10 ns difference in arrival time in the NE213 cells, 
and thus can be treated separately. 

In the 18% of cases where the 7 is detected, most of the time it will leave 
only part of its energy via Compton scattering. The resulting distribution 
of energy deposited in the Ge will be flat, starting at the recoil energy and 
extending up to ~500 keV. 

Multiple Interaction of a Neutron in Ge 

The interaction length of 2.32 MeV neutrons in Ge is approximately 6 cm. 
Despite the small size of the Ge detector chosen for the experiment (3.1 
cm 3 ), incident neutrons have a 20% probability to interact twice or more in 
the crystal. The corresponding coincident rates are approximately half of 
the elastic rate (Table §). It is not reduced by the cut on the time-of-flight 
of the scattered neutron which is not significantly altered by the successive 
interactions. However the angle of the outgoing neutron is uncorrelated with 
the energy deposited in the Ge detector. Some rejection of this background 
can be thus achieved, as single-scattering neutrons detected at a given angle 6 
should correspond to a single E ee value within experimental resolution, while 
the corresponding distribution for multiple scatters is a continuum extending 
approximately from E ee to twice E ee . 

Collisions in the Surrounding Materials 
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Neutrons can also have multiple interactions, with at least one occurring 
in the Ge and the other(s) in the materials surrounding the detector. In the 
SICANE setup these additional interactions occur mainly in the cryostat[| 
and the chamber containing the hydrogenated target. Once the time-of-flight 
requirements are applied, the coincidence rate of such events is predicted to 
be of the order of the elastic rate. However, for half of these events, the energy 
deposited in Ge is only a few keV and will not pass the trigger requirement. 
For the other half, the deposited energy is evenly distributed between and 

E ee 

In order to treat adequately the background of events coming from i) 
inelastic recoils followed by Compton-scattering, ii) multiple neutron scat- 
tering in the Ge itself and ii) multiple scattering in the surrounding material, 
a smooth background is included in the fits to the observed E ee distributions 
yielding the average value at a given angle and time-of-flight selection. 



4 Scintillation Quenching Factor of Nal(Tl) 

As a first test of the performance of the SICANE facility, it was used for 
a measurement the scintillation quenching factor of Nal(Tl). Measurements 
with Nal(Tl) are considerably easier to perform than measurements with Ge. 
There is no need for a cryostat and in addition the good timing resolution 
provides an additional test to remove random coincidences: the neutron time- 
of-flight between the hydrogenated target and the Nal(Tl). 

In order to reduce multiple scattering, the Nal(Tl) detector had a small 
size (25 mm diameter, 25 mm thickness). Its energy response E ee for photons 
was calibrated in the energy range from 60 to 511 keV using 57 Co, 241 Am 
and 22 Na radioactive sources. With incident neutron energy of 2.1 MeV only 
sodium recoils are observed at all four scattering angles of 45°, 90°, 120° and 
165°, 

To simulate the experimental conditions for the Ge experiment, the total 
count rate in the Nal(Tl) detector was limited to 35 Hz, yielding a coinci- 
dence count rate of 0.33 Hz. Figure |5| shows recoil energy spectra recorded 
at different scattering angles after the selections on the neutron identification 
parameter in the NE213 cells and on the total neutron time-of-flight between 
the polyethylene target and the NE213. The largest peak observed in each 

3 The total thickness of the copper cryogenic shields is 4mm. 
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spectrum corresponds to the elastic scattering of neutrons on Na nuclei. The 
peak around 60 keV, clearly seen at #=45°, is due to inelastic scattering via 
the first excited state of 127 I. The structure below 10 keV evolving at the 
different angles could be due to neutron elastic scatterings on 127 I, produc- 
ing a peak expected at approximately 6 keV at #=165°, and at even lower 
energies at other angles. An additional cut on the time-of-flight between the 
Nal(Tl) and NE213 detectors does not remove a significant number of events. 
This indicates that the neutron identification and the total time-of-flight are 
sufficient to select a clean sample of neutron scattering events. 

The En values for Na recoils in elastic scattering for the different scat- 
tering angle values and the corresponding quenching factors values are listed 
in table [I]. As shown in Fig. the measured quenching factors are in good 
agreement with previous results ||. Although the present measurement is 
based on a one-day test experiment, the statistical errors are already below 
the systematic uncertainties. At low recoil energy, the dominant uncertainty 
is the recoil energy due to the finite solid angle of the NE213 cells. At high 
recoil energy, the main uncertainty is the 7 energy calibration. 



5 Ionization and Heat Quenching Factors of 
a Germanium Bolometer 

5.1 Detector principle 

The EDELWEISS collaboration fM searches for WIMP Dark Matter us- 



ing heat-and-ionization cryogenic germanium detectors, also called bolome- 
ters. The energy deposited in the detector produces two signals. First, the 
electron-hole pairs created by ionization are collected on electrodes implanted 
or vaporized on the surface of the crystal. Secondly, a thermistance measures 
the few fiK heating of the crystal due to i) the thermalization of the electron 
or nuclear recoil from the primary interaction, and ii) the Joule energy as- 
sociated with the migration of the electron-holes due to the applied electric 
field. Whether the initial recoil is electronic or nuclear, the subsequent Joule 
heating (or Luke-Neganov effect f27|) depends only on the strength ionization 



signal, and thus is not affected by the quenching factor. Correcting for this 
effect in the presence of a bias V Q on the electrodes, the quenching factors 
associated with the ionization and recoil energy measurements (Q and Q', 
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respectively) are 



Q 



Ej_ 



(5) 



Q 




pair 



(6) 



where V pa i r = 2.96 V is the electron-hole pair creation potential in Ge. Ej 
and Eh are respectively the electron-equivalent ionization and heat signal 
amplitudes, i.e. as calibrated used 7-ray sources. Er is the recoil energy as 
given by Eq. ||. 



For the measurement of its quenching factor, the cryogenic Ge detector is 
placed inside a dilution refrigerator. A base temperature of 25 mK was 
achieved, however it was regulated at 35 mK in order to optimize the prop- 
erties of the heat signal for this measurement. The thicknesses of the different 
screens of the cryostat were chosen to minimize neutron scattering. Special 
care was taken to avoid mechanical vibrations. In order to isolate the cryo- 
stat from the beam line components and the multidetector array elements, 
it is maintained at beam height (1.8 m) with a wooden structure mounted 
on pneumatic insulators. A metallic platform above the cryostat facilitates 
liquid nitrogen and helium fills and access to the bolometer electronics. The 
center of the detector is 60 cm away from the polyethylene target. 

5.3 Detector 

The bolometer used in this work is a cylindrical crystal with a 10 mm height 
and a 20 mm diameter. In this prototype detector [^] the heat sensor is a 
NbSi thin film used as an Anderson insulator. The operating potential V Q is 
6.26 V. To facilitate the detector monitoring, a weak collimated 57 Co source 
is mounted close to it. 

Keeping the data rate as high as possible is desirable in this statistics- 
limited measurement. In this respect, an advantage of NbSi thin film sensors 
relative to the NTD sensors normally used by the EDELWEISS collaboration 
(see e.g. Ref. p6| ) is the relatively fast risetime and decay of the heat signal. 



5.2 Cryostat 
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These time characteristics depend on the temperature of the detector. At 
35 mK, the risetime is approximately 1 ms and the decay of the signal has a 
fast (~4 ms) and a slow (~40 ms) components. The fast component is due to 
ballistic phonons and its amplitude depends not only on the deposited energy 
but also on the interaction point relative to the sensor. The slow component 
corresponds to thermal phonons and depends only on the deposited energy. 
It limits the acceptable interaction rate in the detector to approximately 
15 Hz. 

The electron-equivalent energy resolutions on ionization and heat signals 
in this experiment are 8 and 41 keV FWHM at 122 keV, respectively. 



5.4 Ionization quenching factor of Ge 

The ionization measurement has a trigger threshold of 17 keV. With a quench- 
ing value of approximately 30%, the expected ionization signal associated 
with neutron scattering in the 90° ring is only 18 keV. Consequently, the 
NE213 detectors were grouped on the 120° and 165° rings, corresponding to 
recoil energies of 93 and 123 keV, respectively. 

Fig. shows the time-of-flight spectra between the polyethylene target 
and the neutron detectors. In Fig. [7|a, only those events rejected by the 
n/7 discrimination variable of the NE213 detectors are displayed. The peaks 
associated with photons produced either from 7 or neutrons interacting in 
the Ge detector (labeled 77 and n.7, respectively) are clearly seen. Accepted 
events in the 120° and 165° rings are shown in Fig. [F]b and c, respectively. 
The 77 and ny peak are strongly suppressed. Neutrons associated with 
an elastic interaction in the Ge detector are expected at time-of-flight s of 
77.1 and 77.4 ns, for the two respective rings. The broad structure from 20 
to 40 ns corresponds to neutrons created in the polyethylene target hitting 
directly a NE213 cell without prior interactions. This implies that the energy 
deposit in the Ge detector is not related to the neutron in the NE213 cell, 
and could be related to another particle emitted by another interaction in 
the polyethylene target in the same beam bunch. 

The most important inelastic states expected to be excited are the 563 
keV states in 76 Ge, 596 keV in 74 Ge, and 835 keV in 72 Ge. The corresponding 
time-of-flights are 84.5, 85.1 and 89.7 ns at 120°, respectively, and approxi- 
mately 0.5 ns later at 165°. An excess of counts is indeed observed at these 
time-of-flights in Figs. |7|b and c, especially at 165° where the simulations 
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predict the largest inelastic yields. 

To clearly assess whether these events are indeed due to inelastic scat- 
tering, the events are categorized according to 3 intervals of time-of-flight, 
corresponding to 3 scattering processes: 70 to 80 ns (elastic scattering), 80 to 
87 ns (inelastic scattering with E 7 =563 or 596 keV), and 87 to 92 ns (inelas- 
tic scattering with E 7 =835 keV). The associated ionization spectra for these 
three categories at the two angles are shown in Fig. [| A peak of events at 
low ionization energy is observed in each category. Its mean energy is largest 
in the elastic selection, and decreases with increased E T The peak loca- 
tions are consistent with those calculated (dashed lines) with the expected 
recoil energy for each event category, assuming a common quenching value of 
33.3%. The fact that a flat background is more present in Figs. |8|c to f than in 
Figs. |]a and b provides further support for the identification of inelastic scat- 
tering, as this behavior is expected as emitted photon may deposit some of 
its energy in the Ge detector. However, the relative amount of smooth back- 
ground may also be partly due to the difference in the signal-to-background 
for each event category (as seen in Fig. |7]). 

A final check of the assignment of the events of Figs. [8]c to f to inelastic 
collisions is to look back at the time-of-flight spectra of the neutron-selected 
events, but this time separating the events in two samples, according to 
the recorded ionization signal. Fig. |9|a shows the time-of-flight spectra of 
events with ionization energies clearly associated with the peak observed in 
Figs, ^a and b, that is, between 25 and 50 keV at 120° and 40 to 65 keV at 
165°. Fig. |9|b shows the time-of-flight spectra of events outside this ionization 
energy intervals. As expected, the elastic scattering peak is well suppressed 
in the second sample. 

These tests clearly demonstrate the presence of inelastic scattering events, 
and indicate that the time-of-flight selection is adequate for identifying them, 
given the present statistics. This identification is important since the differ- 
ence in recoil energy between elastic and inelastic scattering events is about 
20%. A improper identification of the elastic sample leading to a contamina- 
tion of inelastic events could have lead to a significant bias in the measure- 
ment of the quenching factor. In our case, assigning the elastic value for the 
recoil energy to all events in the 68 to 92 ns range would have reduced the 
average Q value measurement by 5%. 

It should be noticed that, with sufficient statistics, inelastic events in- 
crease the number of recoil energies covered by the measurements, going 
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down to lower values (see Fig. ||). 



The quenching values extracted from Fig. |] data are listed in table § 
and shown on Fig. |TI]. The statistical precision on the average Q value 
obtained in a one-day experiment is 2%. The results are consistent within 
statistical uncertainty with published results [|ll], [13], |14]] and the Lindhard 
model calculations. No significant systematic shifts are observed between the 
measurements of Ge(Li) in liquid nitrogen or our measurement of hyperpure 
Ge at cryogenic temperatures, or between our direct measurement of Q and 
the Q/Q' measurement of Ref. O] in the same temperature range. 



5.5 Heat quenching factor of Ge 

The analysis for the heat quenching follows closely that of the ionization 
channel. The neutron and time-of-flight selections are identical. 

However, in addition to the significantly worse resolution of the heat signal 
relative to ionization, the slow components of the heat signals are strongly 
affected by pile-up and microphonics. In order to get a heat-to-ionization 
ratio of one for events clearly associated with photons (ionization above 100 
keV and identification in the NE213 cell), it is necessary to impose that 
the energy deduced from the slow and fast components agree to within 75 
keV. This ensures the integrity of the adjusted pulse shape, at the price of 
rejecting all surface events with anomalously large fast components. Clearly, 
the problem of pile-up is a strong limitation for quenching measurements 
with slow cryogenic detectors. 

Given the reduced statistics, a meaningful measurement is only obtained 
for the total sample, dividing each measured energy by the recoil energy 
corresponding to the given angle and time-of-flight selection. The resolution 
does not allow for the identification of those inelastic events where the photon 
was partly detected. In order to be able to use inelastic events despite this 
problem, events with ionization quenching values below 20% and above 45% 
are rejected. 



The resulting distribution of heat quenching values is shown in Fig. |ll|a. 
A Gaussian fit to the distribution gives an average quenching of 50 ± 3 %. 
Subtracting the Luke-Neganov effect corresponding to the applied electric 
field of 6.24 V using our measured average quenching value of 33.3%, we 
obtain (Fig. |TT|p) a quenching factor for the nuclear recoil part of the signal 
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of 87 ± 10%. The quoted error is statistical only; the systematic error, 
mostly due to calibration uncertainties, is of the same order. Within these 
large errors, the result is consistent with the expected value, which should 



be very close to 100% [1|, [17] . 



This one-day experiment clearly illustrates what can be achieved with 
this technique. The present statistical error would scale with the improve- 
ment of resolution in the heat channel. With sufficient energy resolution and 
statistics, a measurement at the percent level may be possible, if systematic 
effects such as inelastic scattering are kept under control. 



6 Conclusion 



The commissioning tests of the SICANE array and its setup confirm its high 
relevance for the calibration of the response of cryogenic detectors to nuclear 
recoils. 

Using inverse reactions to produce a highly collimated neutron beam, no 
massive shielding of the neutron detectors is necessary, giving more freedom 
in their spatial setup. Competitive results have been obtained in a short 
running time (less than a day) for the quenching of scintillation in Nal(Tl) 
and of ionization in Ge. The array also makes possible a deeper investigation 
of the effect of inelastic scattering on the measurement. 

The powerful diagnostic tools provided by the array (simultaneous mea- 
surements at different angles, neutron identification and time-of-flight mea- 
surements) are ideally suited for the study of cryogenic detectors. For such 
measurements, important factors that require proper attention are: i) the 
slow response of the detector, and how it compares with the single rate in 
the cryogenic detector and the rate in coincidence with the neutron array; ii) 
the rate in the neutron detector array due to neutron scattering in the mate- 
rial in the vicinity of the detector and Hi) the ability to resolve inelastic and 
elastic events, either using the energy or the time-of-flight measurements. 

In addition, for the first time, it was directly verified that the ionization 
quenching factor of Ge at 35 mK is very close to that at liquid nitrogen 
temperature. This is consistent with the fact that the CDMS and EDEL- 
WEISS collaborations |T3|, [B5| obtain Q/Q' ratios compatible with the Q 
values measured at liquid nitrogen temperature JTT|, [L4| under the assump- 
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tion that Q'—l. Conversely, this last hypothesis is verified for the first time 
in the present work, albeit at a 15% level. 
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e 




Q 




(keV) 


(%) 


45° 


50 ± 3 


28.1 ± 2.8 


90° 


171 ± 6 


29.4 ± 2.9 


120° 


256 ± 8 


27.5 ± 2.8 


165° 


336 ± 10 


25.4 ± 1.2 



Table 1: Quenching values measured for Na recoils induced by the elastic 
scattering of 2.1 MeV neutrons for the different scattering angle values. 





Coincidence rate 


Process 


(/hour) 






165° 


120° 


90° 


Elastic 


2.0 


1.6 


0.2 


Inelastic 


1.7 


0.6 


0.1 


Multiple scattering inside Ge 


1.0 


0.6 


0.1 


Multiple scattering outside Ge 


1.9 


2.8 


0.3 



Table 2: Simulated rates of Ge-NE213 coincidences after neutron identifi- 
cation and requiring that the neutron time-of-flight from the hydrogenated 
target to the NE213 counter is between 70 and 85 ns. The assumed incident 
yield on the target is 70 neutron/s. 



e 


Process 


E R (Ge) 
(keV) 


Q(%) 
(Stat, error only) 


120° 


Elastic 


94 ± 3 


34.7 ± 0.8 




Inel. (596) 


80 ± 2 


32.8 ± 0.9 




Inel. (835) 


76 ± 2 


30.1 ± 4.0 


165° 


Elastic 


123 ± 4 


33.7 ± 0.7 




Inel. (596) 


104 ± 3 


31.4 ± 1.2 




Inel. (835) 


100 ± 3 


31.4 ± 1.4 




Average 




33.3 ± 0.5 



Table 3: Ionization quenching measured for Ge recoils induced by the elastic 
scattering of 2.32 MeV neutrons, for the different scattering angle values and 
scattering processes. The global systematic error on Q due to calibration 
and beam energy uncertainties is 3% (not included here). 
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Figure 1: Time of flight spectra between the arrival of the n B beam burst on 
the hydrogenated target and the detection of a particle in a NE213 neutron 
detector located at angles varying from 0° to 25°. The incident n B beam 
energy is 35 MeV. 
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Figure 2: Neutron energy as a function of the emission angle for the reaction 
p( 11 B, 11 C)n at 35 MeV beam energy. The line is the expectation from Eq. 
The experimental points corresponds to the peaks observed in the preceding 
figure. The horizontal error bar correspond to the solid angle of the neutron 
detector and the vertical one to the neutron energy spread associated to the 
width of the observed peaks. 
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Figure 3: Time of flight spectra between the arrival of the 7 Li beam burst on 
the hydrogenated target and the detection of a particle in the NE213 neutron 
detector at 6 = 0° for 7 Li beam energies equal to (a) 13.15 (b) 13.4 and (c) 
14.9 MeV. 
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Figure 4: Experimental setup with the four rings (45°, 90°, 120° and 165°) 
of neutron detectors and its support structure. The cryostat and the reaction 
chamber are also shown. 



24 




50 100 150 200 

Equivalent-Electron Energy (keV) 



Figure 5: Energy spectrum (in equivalent-electron keV) recorded in the Nal 
detector in coincidence with a neutron identified in a NE213 cell at 9 = 165°, 
120°, 90° and 45° (a to d, respectively). The number of counts below 24 
keV at 165° has been scaled down by a factor of 3 for clarity. The dashed 
lines show the expected peaks positions for elastic (Na,I) and inelastic (I*) 
scatterings, assuming Q(Na)=0.27 and Q(I)=0.09. 
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Figure 6: Scintillation quenching factor for Na recoils in Nal(Tl) as a func- 
tion of the recoil energy. The experimental data are from Refs. and from 
this work. 
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Figure 7: Distribution of time-of-flights relative to the beam bursts, ob- 
served in the neutron detectors at 120° and 165° when selecting (a) events 
identified as 7, (b) events identified as neutrons at 120° and (c) at 165°. The 
time intervals expected for (7,7'), ( n ,7), ( n , n ) or (n,n') scattering in the Ge 
detector are indicated as dashed lines. Direct neutrons are those going di- 
rectly from the hydrogenated target to the neutron detector without passing 
through the Ge detector. 
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Figure 8: Ionization energy spectra (in equivalent-electron keV) recorded in 
the Ge detector in coincidence with a neutron identified in a NE213 cell at 9 
= 120° (a,c,e at left) and 165° (b,d,f at right). The elastic selection (a and 
b) corresponds to total neutron time-of-flights of 70 to 80 ns, the "inelastic 
(596)" selection, to 80-87 ns (c and d), and the "inelastic (835)" selection, 
to 87-92 ns (e and f). 
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Figure 9: Distribution of time-of-flights relative to the beam bursts, observed 
in the neutron detectors at 120° and 165° for events identified as neutrons, 
with (a) the cut on the ionization signal in Ge that favors elastic events (see 
text) and (b) the cut favoring inelastic events. The expected time-of-flights 
of elastic and inelastic events are marked with arrows. The dashed lines 
represent the three time-of-flight intervals used in the previous figure. 
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Figure 10: Ionization quenching factor measured for Ge recoils in Ge(Li) 
detectors at liquid nitrogen temperature 0, [12|, [Uj, and Ge recoils in hy- 



perpure Ge at 35 mK (this work). The Q/Q' measurement of Ref. [TB| is also 
displayed, assuming a heat quenching Q' of 100%. The thick line represents 
a calculation based on the Lindhard theory [M, as parametrized in Ref. B. 
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Figure 11: Quenching factor for (a) the heat signal and (b) the recoil energy 
deduced from the difference between the heat and ionization signals. The full 
lines correspond to fits to the data yielding the values discussed in the text. 
The dashed lines represent the expected values assuming Q = 33.3% and Q' 
= 100%. 
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